Quantum effects for the 2D soliton in isotropic ferromagnets 
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We evaluate a zero-point quantum correction to a Belavin-Polyakov soliton in an isotropic 2D 
ferromagnet. By revising the scattering problem of quasi-particles by a soliton we show that it leads 
to the Aharonov-Bohm type of scattering, hence the scattering data can not be obtained by the 
Born approximation. We proof that the soliton energy with account of quantum corrections does 
not have a minimum as a function of its radius, which is usually interpreted as a soliton instability. 
On the other hand, we show that long lifetime solitons can exist in ferromagnets due to an additional 
integral of motion, which is absent for the cr-model. 

PACS numbers: 75.10.Hk, 75.30.Ds, 05.45.Yv 



I 

CO 



I 

S3 
O 

o 



(N 
> 

in 

(N 

o 

03 



i 

S3 
O 

o 



X 



Solitons are known to play an important role in several 
branches of field theory and condensed matter physics, 
see Ref. [j] for a review. In particular, solitons treated as 
nonlinear excitations are important in ID and 2D mag- 
netism ||, |j. A serious impediment in studying 2D 
spin systems arises due to the absence of exact analytical 
solutions for most models. Thereupon special attention 
is deserved to models which admit an analytical treat- 
ment. One of the well-known examples is a model of the 
2D isotropic ferromagnet (FM), which provides an exact 
analytical soliton, the so-called Belavin-Polyakov (BP) 
soliton In terms of the normalized magnetization, 
m = (sin 9 cos <j>, sin 9 sin <f), cos 9), the soliton structure 
rriBP is described by the formula || 
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Here p and x are the polar coordinates in the magnet 
plane, the integer q is the ^-topological charge of the 
soliton, R and ipo are arbitrary parameters. BP-solitons 
of the form (|l]) appear in different models of non-linear 
field theory and condensed matter physics . In partic- 
ular, the BP-solitons are important for the ferromagnetic 
quantum Hall effect (|]. 

The unique problem of the BP soliton is that its en- 
ergy Sbp — 4:nJS 2 \q\ (J is the exchange integral, S is the 
atomic spin) does not depend on its radius R: this results 
from the scale invariance of the system, which is a part 
of the general conformal invariance of the model. Since 
the soliton radius is not fixed, one is then free to let it go 
up to the system size, and the thermal excitation of soli- 
tons will break the long-range order S. However, recent 
studies have shown that a quantization of the soliton of 
the classical cr-model, which can be attributed to anti- 
ferromagnets, breaks the static scale invariance @, ||, [§■ 
The natural question whether it works for FMs is still 
open. 

The purpose of our study is to examine the role of 
quantum fluctuations for the soliton properties. We treat 



the problem semiclassically using the one-loop correction 
to the classical soliton energy, originally calculated by 
Dashen et al. |h], [n]] for ID solitons, see also Ref. |lj. 
To generalize these results to the 2D case one needs to 
solve the soliton-magnon scattering problem for 2D mag- 
nets. We show a unique new feature of the 2D soliton- 
magnon interaction, which is absent in ID: the soliton 
acts on magnons not only by some local potential, but 
also in the same way as an effective long-ranged magnetic 
field acts on a charged particle; this essentially changes 
the scattering picture, leading to the Aharonov-Bohm 
(AB) scenario. We state that the AB type of scattering 
is a general consequence of 2D scattering by a topologi- 
cal soliton. We calculate the Casimir energy for the BP- 
soliton in FMs and show that the quantum correction 
can not provide a fixed size for the soliton. Nevertheless, 
we show that long lifetime solitons can exist in FMs due 
to an additional integral of motion, contrary to antifer- 
romagnets. 

The macroscopic dynamics of the classical FM follows 
the Landau-Lifshitz equations 
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sin 6 d t (j) = V 2 6 - sin 6 cos 6 (V<j 
sin d f = -V • (sin 2 6V6) , 



(2) 



where D is the stiffness coefficient of the spin-waves, 
which are characterized by the dispersion law Lu(k) — 
Dk 2 . To analyze the soliton-magnon interaction, we con- 
sider small oscillations of the magnetization m on the 
background of the stationary BP-soliton m,Bp. These 
oscillations can be described in terms of the complex val- 
ued "wave function" * = 9 — 9-qp + isin#Bp(</> — </>bp), 
see Ref. [l3|. For the further analysis it is instructive to 
rewrite the linearized equation for the ^-function in the 
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form of the Schrodinger equation: 



H»£ = —d t V, H = (-»V - A) + V, (3a) 



V = -%sin 2 B p, 
P 



qcos9 BP 
A = e x . (3b) 



The Hamiltonian H has a form which is typical for 
a quantum-mechanical charged particle in the presence 
of a scalar potential V and an additional magnetic field 
with a vector potential A. The appearance of an effective 
magnetic field is a new feature of the 2D soliton-magnon 
interaction, which is always absent in ID systems. Dis- 
cerning this effective magnetic field gives the possibility 
to draw a number of general conclusions about soliton- 
magnon scattering in the 2D case, see below. 

For the system (||) we apply the standard partial wave 
expansion, using the Ansatz: 

^(p,X,t)= ^2 ^m(p) exp (imx - iuj a t + f3 a ) . (4) 

a— (/c,m) 

Here the integer to is the azimuthal quantum number, k 
is the radial wave number, and (3 a is an arbitrary initial 
phase. Each partial wave ip m is an eigenfunction of the 
2D radial Schrodinger equation 
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Here the term linear in m reminds of an effective mag- 
netic field A. The scattering problem can be formulated 
in the usual way. The eigenfunctions for free magnon 
modes have the form ^ J\m\(kp), with an asymp- 
totic behavior ip { m C (kp)" 1 ^ 2 cos (kp — \m\ir/2 — 7r/4) 
when kp 3> \m\: J m is Bessel function. In the pres- 
ence of a soliton the behavior of a magnon solution can 
be analyzed at large distances, p ^S> R. In view of the 
asymptotic behavior U m ~ |m + q\ 2 / p 2 , in the limiting 
case fcp> |m| one has the usual result p3[: 
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The phase shift r\ m contains all information about the 
scattering process. 

The main features of the scattering on a topological 
BP-soliton are caused by the magnetic field A. As an 
analogue of the Zeeman splitting of electron energy terms 
in an external magnetic field, the presence of an effective 
magnetic field breaks the symmetry i] m (k) = rj- m (k). It 
is necessary to take into account separately positive and 
negative m's. 

Sometimes the soliton-magnon scattering problem is 
treated perturbatively using the Born approximation 
pi, which, in principle, can be used for the scattering in 



a magnetic field. However, due to the topological soliton 
properties, A is a long-ranged field 
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when p^s> R, 

which is typical for the AB effect |L4j . For such a type of 
scattering, some standard scattering results fail. For ex- 
ample, the Levinson theorem must be modified for long- 
range potential systems jlB). Since scattering phase 
shifts are not still localized, there appears a problem of 
the regularization of the scattering series like in conven- 
tional AB scattering picture [p^| . As was firstly noted 
by Feinberg [jlTj, the Born approximation fails for the 
AB scattering; it gives an average of two different modes 
with opposite signs of m |ll| . Thus we need a more pre- 
cise approach than the Born approximation. One needs 
to stress that such a long-range behavior is not a result 
of slow algebraic decay of the out-of-plane structure of 
the BP soliton. It is a consequence of the topology of 
the BP-soliton, namely, of the relation <f> = q\, thus the 
AB-type of scattering is valid also for anisotropic mag- 
nets @f§,|2l|. 

Let us discuss the soliton with the topological charge 
q = 1, which has the lowest energy. Such a soliton has 
two internal zero-frequency modes, which are the limit 
of the continuum spectrum as k — > 0: |22 
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The mode with m = +1 is a local translational mode, 
which describes a soliton shift, the mode with to = is 
the half-local rotational mode. The mode with m = +1 
has an exact analytical solution for any finite values of k 
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fpm=i(p) = -h{kp) - 
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kp (p/R) 2 + V 

hence this mode does not scatter at all, ?/ m= +i = 
p3| . Note that in the interesting case of long-wavelength 
asymptotic behavior (kp <C 1) at large distances p ^> R, 
this expression has the same form as a combination of 
Bessel and Neumann functions, J2(kp) oc (kp) 2 and 
Y 2 (kp) oc (kp)~ 2 . Thus, the second term in (g) imi- 
tates the presence of the function Y-x and the presence 
of scattering, which caused a conclusion in Ref. || that 
the mode with m = +1 can be scattered. The scattering 
phase shift can be found in both limiting cases, for small 
and large dimensionless radial wave number >c = kR [ p"3"| . 
For long-wave lengths, x <C 1 
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and in the opposite case of short-wave lengths |T3| 

7r(TO — 1) 

Tj m ~ 7rsgn(m-l) . (9b) 

Now we are able to calculate the density of magnon 
states. Let us generalize the main arguments of Dashen 
et al. (ll], |ll| for the 2D system. The idea of the ap- 
proach is to calculate energy shifts of vacuum magnon 
states in the presence of a soliton, which are constructed 
as one-loop quantum corrections to the soliton energy. 
The energy of the vacuum comes from the zero-point 
fluctuations of the magnon states. Without the soliton, 
each vacuum magnon makes a contribution as /LDfc^ ac /2, 
where {fc V ac} is set of allowable wave vectors. In the soli- 
ton presence the set of allowable wave vectors changes, 
{k}. The energy of the state with the wave vector k is 
hDk 2 /2. Therefore the energy correction is 
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FIG. 1: The change 3? in the density of magnon states due to 
the soliton vs dimensionless wavenumber (k = kR), obtained 
by numerical integration of the Schrodinger equation (|s]). 
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To determine the set of allowable states, we put the sys- 
tem in a very large box of the size L, making all states 
discrete. In the limiting case L — > oo the energy correc- 
tion (|l0|) does not depend on the form of the boundary 
conditions. For the 2D case we choose fixed boundary 
conditions for a circular box of radius L |L3|]. Since the 
free magnons are described by the Bessel function ip m e °' 
by enforcing fixed boundary conditions ip m cc (P = L) = 0, 
we fix the allowed values of the radial wave number, 
k™ c L = jm , where jm** is the n— th zero of the Bessel 
function J m . In the region of interest, n>l, the zeros 

(n) 

of the Bessel function j m » irn. Thus the allowed val- 
ues of the wave numbers are fc™ c « irn/L, similar to the 
ID case. However, the above-used simple equation for 
jljb' is valid only if |m| is not very large. In the case of 
|m| 3> 1, the first zero of the Bessel function j$ w \m\. 
Hence, in a finite system there appears a restriction for 
the allowed number of modes, |m| < L, and the sum rule 
for the 2D case takes the form: 



k.m u m— — kL 



(ii) 



For magnon states in the presence of the soliton there 
appears a phase shift r\ m due to the soliton-magnon scat- 
tering, therefore 



kT c L 



(n) 



V 

L' 



The one-loop correction to the soliton energy reads 
J ,i- 1 oo P = kv (k)dk, v(k)=J2 Vm(k). (12) 



Here rj(k) is a sum of the phase shifts of all partial waves 
with the fixed radial wave number k. Note that our def- 
inition of r)(k), in contrast to the case of the cr-modcl 
contains independent summations over positive and 
negative m, which reflects the breaking of the symmetry 
m — > —m. 

The short-wave length behavior of the phase shift is 
responsible for ultraviolet (UV) singularities. Follow- 
ing Refs. H in order to avoid the UV singularity one 
needs to derive the short-wave length asymptotics for 
the phase shift. There appears a problem to sum to in- 
finity an alternating series, which has no absolute con- 
vergence. Symmetric limits regularize this summation 
similar to the exponential regularization for the original 
AB-scattering |l6| ; different waves are taken into account 

in) 

in the order in which the poles km in the scattering am- 
plitude appear as k increases p4| . Finally, we found that 
77(00) = — 7r for the soliton with q = 1^5|. Note that 
this result can not be obtained from the Born approx- 
imation because of the long-range nature of the AB- 
scattering; namely, perturbative Born calculations re- 
sulted in different conclusions that 77(00) = n (see Ref. |?]) 
and 77(00) = 2tt (see Ref. ^). 

To check our analytical predictions we have calculated 
77(fc) by the numerical integration of the Schrodinger 
equation (||) for rn 6 [—100; 100], which gives 77(00) = — tt 
with the precision 10 -3 . 



Finally, the one-loop correction reads S 1 loop = S( 
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Here the counterterm <?ct = — ^■'n(k)k 2 1~ can 
be evaluated using a momentum cutoff A, which results 
in Jot = ^HDA 2 . Taking away this UV term, e.g. by 



renormalizing the exchange constant 0], we ends with 



4 



the finite Casimir energy: 
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Here describes the change in the density of magnon 
states due to the soliton. This expression can be ana- 
lyzed analytically in limiting cases. Using the asymp- 
totical behavior for the phase shift for different modes 
one can conclude that the maximum scattering in the 
long wavelength limit corresponds to the mode with 
to = 0, hence the density of states has a singularity: 
!R(x) ~ — (2x) _1 In 2 x — > — oo when k — > 0. In the short 
wavelength limit all modes compensate each other and 
— > oo) — > 0. In the intermediate range the change 
in the density of states can be found numerically only, 
see Fig. ^|. One can conclude that for all >c the change 
in the density of states takes only negative values and 
J D?(>f)d>f = —1. Finally, the Casimir energy is ob- 
tained as 
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C= hS\Jl(x)\dx, (14) 



where C is a constant, which we calculated numerically 
to be C ~ 0.38. We introduced in ( p"4| ) the typical length 
scale a = \JhDj JS, which is about a lattice constant. 
Note that this parameter is absent for the static BP- 
soliton problem, but it naturally breaks an initial scale 
invariance of the model in the dynamics, see Eq. (||). 
The soliton energy is reduced when the soliton radius 
decreases. 

Let us discuss physical consequences of Eq. ([lj]). First, 
the soliton energy with account of the quantum correc- 
tion does not have a minimum as a function of the soliton 
radius. Usually this is interpreted as a soliton instability 
in the context of the Hobart-Derrick theorem, see Ref. |. 
We will show here that this property leads to a dissipa- 
tion of the soliton energy caused by magnon radiation, 
common to that for 3D Hopf solitons in isotropic FMs 
p3| . As an important contrast to cr-models, a model 
of the FM has an additional integral of motion, the z- 
component of the total spin S z ~ 2irS{R/ a) 2 [Q. Even 
for the static limit S z takes a nonzero value; the energy 
dissipation caused by the radiation of magnon pairs with 
wave vectors k and — k will be accompanied by a decrease 
of the value of S z by two. Therefore the soliton lifetime 
r = S z /(dS z /dt) can be sizeable, when S z is large (S z is 
the number of bound magnons in the soliton). 

The amplitude of the radiation process is w ~ 
JS(ak) 2 



see 

dS. 
df 



In accordance to Fermi's golden rule, 
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Here E is the total energy of the soliton with account of 
the Casimir energy. Calculating dS z / dt in the continuum 



limit, one can finally obtain the soliton lifetime 

AttH [R\ W , N 

T ~J(?(a) ■ (15) 

Note that the lifetime is much bigger than h/J for R > a. 

To conclude, quantum effects decrease the energy of 
the BP soliton in isotropic FMs, more strongly for small 
soliton radius. Nevertheless, the original argumentation 
by Belavin and Polyakov || about the breaking of the 
long-range order of the system is still valid. It is based 
on the fact that in isotropic FMs the energy is indepen- 
dent of the soliton radius R, and R can be comparable 
with the system size, R ~ L. However, in the case of 
large radii the quantum correction is negligible, the en- 
ergy has a finite limit for R — > oo, so the problem of 
long-range ordering has the classical form. Another as- 
pect of the problem is the fate of the BP soliton at small 
finite R. As can be seen form Eq. (|l5|), the lifetime is 
small for small R. In some respects the situation is sim- 
ilar to the problem of the black hole evaporation, i.e. 
the large radius soliton dissipation is very slow, and can 
be neglected. At the same time the dissipation of small 
radius solitons is very fast. When the soliton radius is 
small enough, the speed of dissipation increases rapidly; 
in the final stage with S z ~ S the soliton can collapse by 
a quantum jump, which is accompanied by a change of 
the topological charge. 
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